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Abstract This report presents a numerical procedure for estimating the true scattering 
coefficient, <r°, from measurements made using wide-beam antennas. The use of wide-beam 
antennas results in an inaccurate estimate of <r° if the narrow-beam approximation is used in the 
retrieval process for cr°. To reduce this error, we propose a correction procedure that estimates 
the error resulted from the narrow-beam approximation and uses the error to obtain a more 
accurate estimate of <r°. An exponential model has been assumed to take into account the 
variation of cr° with incidence angles, and the model parameters are estimated from measured 
data. Based on the model and knowledge of the antenna pattern, the procedure calculates the 
error due to the narrow-beam approximation. The procedure is shown to provide a significant 
improvement in estimation of cr° obtained with wide-beam antennas. The proposed procedure is 
also shown insensitive to the assumed <r° model . 
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1. Introduction 


Scatterometers, calibrated radar systems, are designed to estimate, accurately and precisely, the 
backscattering coefficient <j° of homogeneous targets from the received power. Because the 
backscattered signal is the superposition of all scatterers illuminated by the antenna beam, the 
implicit assumptions of this retrieval procedure are that those illuminated scatterers have the 
same property as far as backscattering is concerned, and the radar system has the same response 
to all the backscattered components. Practically, these assumptions are approximately valid only 
when a narrow-beam antenna is used in the measurement, and common retrieval procedures are 
based on so-called narrow-beam approximation, which will be explained in following discussion. 


For a monostatic radar, the received power from a distributed target is, mathematically, related to 
<r° by an integral equation, the radar equation, as expressed by eq. (1) 


PrW 


Pr£Gl r 
(4/r) 3 J 


R\e a ,<Pa) 


a\e a , 4 > a )dA 


( 1 ) 


where 

P, = transmitted power 

A = wavelength 

G 0 = maximum gain of antenna 

g 2 = normalized two-way antenna gain function 

o° = scattering coefficient of the target 

R = range between antenna and target 

6 — incidence angle of EM waves (see Fig. 1), 

The radar illumination geometry is illustrated in Fig. 1. Notice that o , g , and R are the 
functions of radar coordinates ( 0 a , <p a ), and the integration is carried out over the illuminated area 
A. The received power P r is a function of incidence angle 0 , as indicated in (1). 


Our objective is to retrieve <r°(0) from P r {6) based on (1), and assumptions are made to facilitate 
the retrieval procedure. One of these assumptions is the narrow-beam approximation, which 
states that, if the antenna beam is sufficiendy narrow, <r°, g 2 , and R can be considered 
approximately as constants over the illuminated area A. Under the narrow-beam approximation, 
variables <r°, g 2 , and R in (1) are replaced by their representative values evaluated at incidence 
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angle 0, and (1) can be simplified as 


P r {0) 


/>,A 2 Gfc°(0) 

(4;r) 3 /? 4 (<?) 


Ai u(0) 


( 2 ) 


where A iU (0) = J dA is the "effectively" illuminated area. A common practice is to define 

A 

Amid) as the area of the effective antenna beamwidth intercepted by the target plane. For a pencil 
beam, the intercepted area is an ellipse. Using (2), cr° can be retrieved from measured data, radar 
parameters and geometry, 


tr?(0) = 


/> r (g)(4jr) 3 *V) 

P,tfGlA in {d) 


(3) 


Equation(3) gives so-called measured <r° (denoted as <r»), which is calculated under the narrow- 
beam approximation. Using (1), (2) and (3), cr?(0) can be related to true <r°(0) as 


<r?(0) = a°(0) 


1 


Aui(0) 
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J R*(d a ,K,d) 
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where 


R{e a ,<t> a ,0) 


R^Pa) 

R(0) ’ 


normalized range 


<r°(0 a ,* a ,0) 


<r°W ’ 


normalized scattering coefficient 


g 2 (0 a , Pa* #), nonnalized two-way antenna pattern as in (1). 


From (4) it can be seen that the cr* represents only an averaged value of the true <x° weighted by 
the normalized antenna pattern and range. However, the condition of the narrow-beam 
approximation ensures that the integrand of (4) is approximately equal to one over the integration 
area, and cr?, the weighted average of cr°, is close to the true cr° at incidence angle 0. 

Unfortunately, wide-beam 1 antennas are often used in scatterometers for various reasons, and the 
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narrow-beam approximation is then not a valid assumption. If eq. (3) is still used in the retrieval 
process, a* could deviate from the true cr° by a large value. Ulaby et al. [1, pp. 755] briefly 
discusses the effect of assuming constant cr° in the illumination integral. With a beamwidth of 
10°, the narrow-beam approximation could result in an error of several dBs in estimating <r° if 
true <r° changes rapidly over the illuminated area. 

Two examples of the effects of the narrow-beam approximation are shown in Fig.2. The 
theoretical values of cr° of a smooth sea-ice surface at several angles are taken from [2, pp. 
1754], and a smoothly fitted curve is generated from those data and extended to the entire 
relevant range of incidence angles. The curve is used to calculate received power from (1) with a 
pencil-beam antenna of Gaussian pattern. The calculated power is then used to calculate a* from 
(3), which is derived under the narrow-beam approximation. As a quantitative illustration, the 
figure shows how large the narrow-beam approximation error could be if a wide-beam antenna is 
used to measure <r° when it changes rapidly over the angular extent of the antenna beamwidth. 
Because the cr° of targets with a smooth surface often exercises a sharp change for small 
incidence angles, it is necessary to reduce the error to an acceptable extent if the narrow-beam 
approximation condition is violated in the measurement. In other words, some correction is 
needed on a? to take into account the effect of average in (4), so that the requisite accuracy of 
estimation is achieved even though wide-beam antennas are used. 

In section 2, we discuss the possible solutions and propose our approach. In sections 3-5, we 
describe the proposed procedure in detail. In section 6, we present two examples, and in section 
7, we give a summary along with our conclusion. 


1. In this report, an antenna is considered as "wide-beam" if the error in the estimate of <7° by eq. (3) is too large to 
be accepted. Otherwise, it is narrow-beam. These terms really involve three factors: antenna beamwidth, 
variation of <r°, and error tolerance in estimating a 0 . 
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Fig. 2 Effects °f narrow-beam approximation: smooth surface of sea ice 



2 . Correction of the Narrow-beam Approximation Error 


Given received power ( P r (6 k )) measured at incidence angles {0*}, we try to estimate true 
<r°(0 k ) even though the narrow-beam approximation may be violated. The difficulty of the 
problem is that, without the narrow-beam approximation, it involves an integral equation (see 
(1)). Observing (1), we can consider the function g 2 (0 a ,^a^B)/ R*{0 a ,^ a ) in integrand as a 
two-dimensional spatial filter. The two-dimensional convolution, as represented by (1), states 
that the input to the filter is the true <r°, and the output of the filter is the power received by the 
antenna (notice that output is one-dimensional, or can be thought as having rotational symmetry). 
We want to find the input signal of the filter from the output, a deconvolution problem. 

For a general two-dimensional deconvolution problem, an inverse filter can be structured to 
retrieve the input signal based on output data as well as characteristics of the filter and the 
measurement noise [3, pp. 61]. However, in our case the facts that the measurements are 
generally made over a small set of incidence angles and that the filter is space variant make the 
use of inverse filters questionable. 

Axline [4] developed a matrix inversion procedure to invert the integral equation. One drawback 
of this method is that it needs P r (0 k ) to be measured with small angular increments. If a 
theoretical or empirical model can be used to describe a 0 , it is possible to develop algorithms to 
estimate the parameters of the assumed model [5] [6]. 

Another approach for attacking the problem is, instead of solving the integral equation (1) 
directly, trying to improve cr,, the estimate based on the narrow-beam approximation, by adding 
a correction term, in terms of dB, which represents the error due to the narrow-beam 
approximation in the retrieval process. In the following discussion, we call the error term a 
narrow-beam approximation error, and name this kind of approach a narrow-beam 
approximation plus correction. The narrow-beam approximation error can be mathematically 
defined as 


£«?*) = 7 .°( 0 *) 


( 5 ) 


A more accurate estimate 


of cr° can be obtained by adding this error term to cr* , 
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t°(0 k ) = <T°(0 k ) + £(0 k ) 

where £(0 k ) is the estimate of E(0 k ) and all terms in (5) and (6) are in dB. 


( 6 ) 


E(0 k ) is unknown and eq. (5) cannot be used to calculate E(0 k ) because <r°(0 k ) is unknown. 
However we can assume that true <?° belongs to (or is close to) a certain model and the model 
parameters can be estimated from measured data. Once the model is determined, the narrow- 
beam approximation error could be calculated via (4). The idea behind this approach is that what 
we need is the error £(0 k ), the difference between &*(0 k ) and a°(0 k ) in (4) due to the effect of 
wide antenna beam, not the model itself, so an accurate model is not crucial as long as the error 
can be estimated from the model with accepted accuracy. In this report, we propose a numerical 
procedure to estimate the narrow-beam approximation error. Later we will see that the proposed 
procedure is indeed insensitive to assumed model. 

The approach of narrow-beam approximation plus correction has been used by several 
researchers, e.g., Kim et al. [6] and Stiles et al. [7]. Kim et al. proposed an exponential model 
and a recursive algorithm to estimate the model parameters. Stiles et al. described an empirical 
formula, based on a group of hypothetical theoretical cr° curves, to calculate the error E(9 k ) for 
small incidence angles. The proposed approach in this report is more general and systematic. 


3. Exponential Model: A Piecewise-Linear Approximation 


In order to estimate the narrow-beam approximation error via (4), a parameterized <r° model as 
the function of incidence angle must be chosen. The choice of the model is based on the 
following considerations: 

a. It fits experimental data, at least approximately. 

b. It facilitates mathematical handling. 

c. It is robust; that is, deviation of true cr° from assumed model should not cause 
unacceptable error. 

Note that the purpose of introducing the model is to estimate the error resulting from the narrow- 
beam approximation, so it is possible to obtain a good estimate of the error even though the 
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model itself may not accurately describe the true scattering characteristics of the target. 2 
A reasonable model is the piecewise-exponential model 


a°(0) = A i txp(-e/B i) e Li _ x ze<e Li , / = 1 n (7) 

Generally N = 2, and the interval &l,) can be identified by inspecting a* or "optimum 
fitting." 

There are two reasons for choosing the exponential model. First, the exponential model has been 
found to fit reasonably well experimental data from North American agricultural terrain, sea ice, 
and ocean[l, pp. 577]. Second, the model is simple and the model parameters are easily 
estimated from measured data. Particularly, the two model parameters, A, and B h are separable 
so that a combined approach of a one-dimensional lookup table plus one-dimensional 
optimization can be established to estimate the parameters, as proposed in this report. 
Mathematically, the exponential model is nothing but a piecewise-linear approximation of the 
true < 7 ° in the log scale (dB). Therefore, we can expect a better model by dividing the covered 
range into more segments, and each segment is modeled with distinct exponential parameters. 
As mentioned earlier, what we are concerned with is the narrow-beam approximation error rather 
than the model itself, and we will give an example to show that even if the true < 7 ° deviates from 
the model, the results of the procedure are still reasonably good. 

As indicated in (7), the so-called exponential model here is actually a piecewise-exponential 
model, or piecewise-liner model in log scale; that is, the model consists of several segments and 
each segment is modeled by a distinct exponential function. For most measured data, at least two 
segments are needed, one covering small incidence angles and another for large angles. Here the 
model is in contrast to the one used by Kim et al., who used one exponential function to cover all 
the relevant incidence angles. 


2. The narrow-beam approximation error is defined in eq. (5), and the estimated error from the model is given by 
eq. (18). Any modeling error will result in errors in cr^, c (d k ) (estimate of <t°( 0*)4) and in <r°(0*) (estimate of 
a?(0 k )). However if the errors are positive correlated, they may cancel each other in the estimate of the narrow- 
beam approximation error given by eq. (18) 
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4. Estimation of Model Parameter : Lookup Table Method 


4.1 Calculated a 0 


With the chosen model, denoted as cr?, the so-called calculated cr® (denoted as cr?) can be 
derived as (see (4)) 


<#**) = 


gV*) f 

Aui(0k) * 


R*(0a,t a ) 


<r 0 „(e a j a )dA 


( 8 ) 


where A iu (S k ) is the equivalent area as in (2). 

As indicated by its name, cr? is obtained by numerically solving eq. (8). On the other hand, cr? is 
calculated from measured P r based on eq. (3). The integration in eq. (4), which relates a* to true 
cr®, is really embedded in the physical process of the measurement although (4) and (8) have 
exacdy the same form. Another difference between (4) and (8) is that cr 0 in (4) is associated 
with the measured target and is unknown, and cr? is the parameterized model that we are trying 
to fit to measured data in the hope that narrow -beam approximation errors regarding cr? and cr? 
are close to each other. It is very clear that all other quantities appearing in (8), including 
antenna pattern and range, should represent the same quantities involved in the measurement 
process. Though errors are inevitable in modeling these quantities, we ignore these errors in 
following discussions. In other words, values of cr* are experimental data; and cr?, on the other 
hand, is only a mathematical expression that reflects the effect of the illumination integral on cr® . 

To avoid confusion, some words are needed here about <r°, cr®, cr?, and cr?. cr® is the true 
backscattering coefficient, which is never known and which we are trying to estimate, a* is the 
measured backscattering coefficient based on received power and the narrow-beam 
approximation (see (3)); it actually includes all effects of noise and errors in the process of 
measurement and calibration, cr? is the estimate of cr® based on the assumed model and 
estimated parameters. The last one, cr?, the calculated backscattering coefficient, is given by (8), 
and it reflects the illumination integral process in (1), which is ignored by the narrow-beam 
approximation. 
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One interesting point is that, even though we are trying to find the best fit of measured data 
within the chosen model, we do not attempt to use < 7 ° as a final estimate of cr°, but only use it to 
find the narrow-beam approximation error. This is because we are not sure whether or not the 
model (and/or the estimated parameters) is accurate. However we believe that the error estimate 
obtained using the model is reasonably accurate for applying the correction procedure. Fig. 3 
gives a concise illustration of the above discussion. Next, we discuss the parameter estimation 
procedure in detail. 


4.2 Problem description 


Comparing (4) with (8), we can see that cr° will equal cr® if everything is a perfect match, 
particularly between the parameterized model cr® in (8) and true <7°, and all other error sources 
are ignored. Therefore, the difference between them can be used to assess the correctness of the 
parameter estimation. The problem can be mathematically described as follows: 

Given a measured vector <r* = (cr?(0|), • • ■ , o? ($„))', where superscript t means the transposed 
vector, over a set of incidence angles { 0j } , we need to find the optimum model parameters A 0 
and B 0 that minimize the squared Euclidean distance between the two vectors, <r? and cr®, 
defined by 


D 1 = ll<7? - a® II 2 (9) 

= £(*?(*) -<£(*,))* 


where <r® is the so-called calculated vector 

whose components { < 7 ®( 0 t ) ) are obtained from (8) with 
<7®(0) = i4exp(-0/fl) 

In (9), the error is calculated in dB because <7° changes drastically with incidence angles. In 
addition, some kind of weighting can be incorporated into the expression to emphasize part of 
the measured data. 
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Fig. 3 Block diagram of correction procedure 




4J Lookup table method 


Various two-dimensional minimum-finding procedures can be used to find the minimum solution 
of (9), but three factors must be taken into account: 

a. Does the procedure converge each time? 

b. Does it converge to the global minimum each time? 

c. How fast does it converge? 

It is not always easy to answer these questions. One way to get around the above difficulties is 
the lookup table method. That is, calculate [<r+(6j)) over possible ranges of parameters A and B 
and store the results in a table in advance. Then in the retrieval process, the optimum parameters 
of the model can be obtained by simply comparing measured data o* with each entry of the table 
and finding the one that has the smallest distance. Obviously the lookup table method can only 
yield a sub-optimum solution in the sense that the parameters A and B must be quantized into 
finite sets so that the minimum solution is "optimum" only as related to the chosen parameter 
sets. The quantized error can be reduced by using a finer quantization step, or a larger-size table. 
On the other hand, we have assumed that the proposed approach is not model-sensitive and we 
are not claiming the exponential model always accurately represents the true <7°. That means the 
narrow-beam approximation error cannot be very sensitive to the model parameters either, so 
that the parameters A and B do not need to be quantized with very small steps. Otherwise, the 
whole procedure is questionable. The incidence angle 9 must also be quantized into a finite set, 
which is less problematic because measurements are generally taken over a set of pre-designed 
incidence angles. Furthermore, interpolation techniques can be used if occasionally the measured 
data have different angles from those in the table. Another problem is the memory requirement 
for storing the table. Totally there are N, x N A x N 0 values in the table where N 0 , N A , and N B 
arc the sizes of quantized 9, A, and B respectively. This requirement may impose serious 
problems in practice. 

One reason for us to choose the exponential model is its mathematical advantage. In particular, 
the two parameters of the model are separable, and the "best" parameter A is easy to obtain if B 
is fixed. As a result, instead of a two-dimensional table, only a one-dimensional table with B as 
parameter is needed for our model fitting, as seen soon. 
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Observe that in equation (8), parameter A of the exponential model only acts as a scale factor. 
Therefore, we can separate the effects of A and B by expressing <r?(0*) as A crt{6 k ) where er_(0*) 
is still calculated from (8) but instead of using the model in (7), the normalized model 


<r°(0) = expH?/B) 
is used in the expression. 

With the above notations, the problem can be restated as 
min D 2 = minllcr? — <r+ll 2 

A,B A.B 

ss min ( min Her? - <r+ll 2 ) 

B A 

- min ( min Her? - Aer?ll 2 ) 

B A 

The significance of (11) is that the two-dimensional search has been translated to two one- 
dimensional searches. This makes it suitable to use only a one-dimensional lookup table to store 
<r°, which is only the function of the model parameter B, and to estimate A directly in the model- 
fitting process. Therefore we not only drastically reduce the storage requirement of the lookup 
table method by a factor N A , the size of quantized A, but also eliminate the quantization error of 
A because now it is estimated over a continuous base. 


4.4 Generation of lookup table 

Each element of the lookup table consists of a particular value of B, say B t out of the set of 
quantized ( B k } , and two corresponding vectors, cr°, and 

The vector <r° ; is defined as 

,.(<£,(*, (12) 

where a%(0 k ) is calculated from (8) using the normalized model in (10) with parameter B r 
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£, is the error vector whose the k lh element, e,(0 k ), is simply the difference between <J°,(0 k ) and 
exp (-d k /Bi), the value of the model at 6 k , both in dB. This error vector will be used for the error 
correction. 

How to choose the quantized parameters set {£, , i = 1,2,...,M} depends on the allowable 
memory size, the acceptable time spent on the calculation, and the tolerable error due to 
quantization on B. The quantization errors can be assessed by observing the difference between 
the components of two vectors, <7.?, and <r° /+1 , corresponding to two adjacent quantized values, B, 
and Bj + The set of [0 it /» 1, . n) is chosen to include all possible incidence angles, in 
consistence with the designed experiment. 

For each quantized B,, <t°,( 0*) is found over the chosen incidence angles {0j\ by numerically 
solving the integral (8). Besides the normalized exponential model (10) with specified B,-, the 
normalized two-way antenna pattern g 2 (. ) must be specified in order to carry out the integral, 
which could be a fitted function of measured antenna pattern or just measured antenna data if the 
measurements were made over fine-enough angles. For the latter case, interpolation might be 
needed in the calculation. The chosen integral area A should be large enough so that the 
contribution to the integral from the outside of the area can be ignored. Note that only the 
relative distance R(0 a , <p a )/R(0 k ) appears in (8), so the real height of the antenna position is 
irrelevant and any height can be assumed for the calculation. 

A Pascal program for generating a required lookup table is listed in the Appendix. The 
integration program used was taken from [8] and slightly modified for two-dimensional 
integration. In the program, the two-way antenna gain was assumed to be a Gaussian pattern. 

Generating the lookup table is a time-consumed task. As a matter of fact, the integration (8) 
needs to be carried out N 9 x N B times where N e and N n are, respectively, the sizes of quantized 
0 and B. However, the table has been made once and for all as long as the same antenna is used. 
The model fitting costs little time once the table has been made. 


4.5 Parameter estimation 


The parameter estimation procedure is straightforward once the lookup table is made. Fig .4 
diagrammatically illustrates the procedure. The squared distance between a* and <7+ is plotted in 
the form of equi-distance contours with parameters A and B as independent variables. Parameter 
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B has been quantized into a set of discrete values { B h i= and a lookup table has been 

made whose ith element is a vector < 7 °, given in (12). For each discrete B, (equivalently, a° ( , 
the element of the lookup table), the local minimum distance Df is calculated with the fixed B, 
and the free parameter A. By comparing [D ] }, the model parameters are declared as the ones 
that yield the smallest distance D* among [Df). As seen from the figure, quantization on 
parameter B results in some error in locating the global minimum. The error can be made as 
small as we wish by reducing the quantization step with the price of a larger table and longer 
searching time. Also two-dimensional interpolation techniques can be used to reduce the error 
further. The step-by-step procedures for estimating the model parameters are as follows: 


Step 1 Calculate <rS vector from measured data using (3). 

Step 2 For each element of the table, <r°, vector, find optimum A h which minimizes the 
distance between the measured vector <r? and cr®^. = according to the 

following equation 

= t £(*?*;)-<<,(*/)) (dB) (13) 

A. jm 1 

which can be proven easily, as shown in following analysis. 

For a fixed B,-, we want to find the optimum parameter A; that minimizes the 
squared distance 

Dj = Her? - A, <7®, II 2 (14) 

= min Her* - Aat ' ll 2 
a 4 

Notice that in the above equation, <r°s are in dB as mentioned earlier, i.e., 

(A<)^ = (A)^+(a®,U 

SO 

Dj — min II <r.°(dB) - A(dB) - a^(dB) II 2 (15) 
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= min t{d(0j)-A-<rlA9j)) 2 W) 

A ;=1 

Differentiate (15) and let the result equal zero, we get (13). 
Step 3 Compare { D ] } to find the smallest one 

D* - min { D ]} 


( 16 ) 


and the corresponding parameters A and B , which yield D , are the fitted 
parameters of the measured vector cr* based on quantized parameters {5;}. The 
error vector corresponding to B can be read out from the table for error correction. 


5 . Estimation of Backscattering Coefficient o° 


Generally, the angle response of cx* can be divided into two regions (see (7)), by either inspecting 
{ a O( 0 .)j or trying different segmentations and finding a suitable division that yields minimum 
combined D*. After the model parameters for each region are determined, the corresponding 
error vectors can be used to correct the narrow-beam approximation errors on the measured data 
{<x?(0, )}. However, the correction at those incidence angles near the intersection of the two fitted 
model segments might not be as accurate as for other angles. A better procedure is to use the 
combined model to directly calculate the error vector. With the fitted parameters, the combined 

model is 



A* exp (~$/B\) 
A\exp(-d/ B\) 


Q <> 6 < do 
do <>e <nfl 


(17) 


The combined model is again used in ( 8 ) to calculate the error of the narrow-beam 
approximation, E(d k ) 
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E(d k ) = <7° mc (e k )-<7° + (e k ) ( dB ) 


(18) 


where <r+(O k ) comes from (8) using the combined model (17). Because E($ k ) reflects the effect 
of the illumination integral (1), we can expect a better estimate of <r° to be 


d 0 (8 k ) = a 0 .(d k ) + E(8 k ) (dB) (19) 

We call d° as corrected <r° because unlike < 7 ? (uncorrected), a correction factor, E(0), has been 
added to overcome the error of the narrow -beam approximation. 


6. Results 


The proposed procedure has been tested on different types of data, and in this section two 
examples are given to show the performance of the procedure. In the first example, the procedure 
was used on data obtained by sampling the curves labeled "measured <r°” in Fig.2 from 0° to 50° 
with steps of 2.5°. The corrected cr 0 s ({<?V*)}), the measured <r°s ({<r?(0*)}), as well as the true 
<r°s, are plotted in Fig.5. The figure clearly shows the improvement in the estimation of <7° using 
the procedure described. Notice that the two-segment exponential model is not a suitable model 
for the smooth surface < 7 ° curve, as shown in the Fig.2. The remarkable improvement justifies 
the earlier statement that the procedure is not model sensitive. 

Another example uses true experimental data that were collected by two different systems. 3 The 
measured <7° (cr*) of both systems and the corrected a 0 (d°) for the wide-beam antenna data are 
plotted in Fig.6. No correction was made on the narrow-beam antenna data because, except at 


3. Two radar systems were used to record data simultaneously during the CRREL 1989 experiment. One radar was 
a FM-CW system with an antenna of 3° beamwidth, and another was a step-frequency system with an antenna 
of 15° beamwidth. 
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Fiffi, 5 Performance of correction: smooth surface of sea ice 




15b 


Fig. 6 Comparison of nanuw and wide beam data (C band, 1-17-89) 



normal incidence, the narrow-beam approximation errors were very small and negligible. The 
figure shows that, as expected, the wide-beam antenna system only gives an erroneous estimate 
of cr° for small incidence angles if the narrow-beam approximation is used to retrieve <r°. On the 
other hand, data from both systems show good consistency after the correction was made on the 
wide-beam data. 


7 . Summary 


In this report, we have presented a correction procedure for the estimation of <7° for data 
collected by a wide-beam system. The procedure estimates the error incurred by the narrow- 
beam approximation and provides a more accurate estimate of <r° by correcting the error in some 
degree. An exponential model is used in the procedure, and a combined approach of a one- 
dimensional lookup table plus one-dimensional optimization is proposed for model fitting. The 
procedure is shown to provide significant improvement on the estimate of <r° and to be 
insensitive to the model assumption of a 0 . 
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9. Appendix Pascal Program For Generating Lookup Table 


PROGRAM BEAMINT; 

CONST 

PI = 3.1415926536; 

HEIGHT =1; 

H2= 1; ( H2 = H*H ) 

D_R = 0.0174532925; 

BEAMWBDTH_3DB = 0.2617994 ; { 15 DEGREES 5 GHZ } 

{ BEAMWIDTH_3DB = 0.0994838 ; 5.7 DEGREES 9.6 GHZ } 

{ BEAMWIDTH_3DB = 0.0877028 ; 5.025 DEGREES 13.6 GHZ } 

S_THETA = 2.5; (DEGREE) 

TYPE 

GLN ARRAY = ARRAY[1..10] OF DOUBLE; 

VAR 

THETA,HALFBEAM_INT,CENTER_YLONG_AXIS,SHORT_AXIS,X1,Y1,Y2,YO : DOUBLE; 
RHO 1 JRH02 J*HI 1 ,SINT,S 1_CIRCLE,S2_CIRCLE^,R0,X_ Y : DOUBLE; 
R_THETAJ > AR ) S_PAR^4.SIGMA r HALFBEAM_EFF > AILL : DOUBLE; 
dRCLE,II,JJ ; INTEGER; 

FDAT : TEXT; 

FUNCTION DB(X:DOUBLE):DOUBLE; 

BEGIN 

DB := 10*LN(X)/LN(10) 

END; 

FUNCTION TAN(X:DOUBLE):DOUBLE; 

VAR 

C: DOUBLE; 

BEGIN 
C CO S(X); 

IF (C * 0) THEN TAN := 1E36 
ELSE TAN ;= SIN(X)/C 

END; 


( arccos function ) 

FUNCTION ACOS(C;DOUBLE):DOUBLE; 
VAR 

SAC : DOUBLE; 



BEGIN 

IF (ABS(C)>1) THEN BEGIN 
IF (C>0) THEN C:= 1 
ELSE C:=- 1 

END; 

IF(C=0) THEN 
AC := PI/2 
ELSE 
BEGIN 

S := SQRT( 1 .0-SQR(C)); 

AC ;= ARCTAN(S/C); 

IF (AC < 0) THEN AC :» PI+AC 
END; 

ACOS :« AC 
END; 


{ antenna gain function ) 

FUNCTION ANTENN A_PATTERN(PHI : DOUBLE):DOUBLE; 
CONST 
C= 1.6651092; 

BEGIN 

ANTENNA_PATTERN := EXP(-SQR(C*PHI/BEAMWIDTH_3DB)) 
END; 


( sigmaO function : exponential model } 

FUNCTION SIGMA0(PHI,PAR : DOUBLE):DOUBLE ; 
BEGIN 

SIGMAO := EXP(-PHI/PAR) 

END; 


{ find angle between the antenna Iook_diiection and ground a pixel } 

FUNCTION ANGLE(RHOJ > HI f R 1 YO,RO : DOUBLE):DOUBLE; 
VAR 

CPHI :DOUBLE; 


BEGIN 

CPHI := (RHO*COS(PHD*YO+H2)/(R*RO); 
ANGLE := ACOS(CPHI) 

END; 


( ANTENNA GAIN function :g**2 } 

FUNCTION FUNC(Pm,RHO:DOUBLE):DOUBLE; 
VAR 
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ANG :DOUBLE; 

BEGIN 

ANG := ANGLE(RHO,PHI,R,YO,RO); 
FUNC :=ANTENNA_PATTERN(ANG); 
END; 


{ Bnd parameters of illumination ellipse } 

PROCEDURE ELUPSE(THETA,HALFBEAM_INT : DOUBLE ; 

VAR CENTER_YLONG_AXIS ,SHORT_AXIS ; DOUBLE); 

VAR 

A,B,C,D,E : DOUBLE; 

BEGIN 

A := SQR(TAN(HALFBEAM_INT)); 

B := SQR(COS(THETA))-A*SQR(SIN(THETA)); 

C := HHGHT*(1+A)*SIN(2*THETA); 

D := SQR(HEIGHT)*(A*SQR(COS(THETA))-SQR(SIN(THETA))); 

CENTER_Y := C/(2*B); 

E := D+SQR(Q/(4*B); 

LONG.AXIS := SQRT(E/B); 

SHORT_AXIS := SQRT(E) 

END; 


( and angle of a ellipse given range in polar coordinate system } 

FUNCTION PHI_LIMIT(CENTER_Y, LONG. AXIS, SHORT_AXIS ( RHO : DOUBLE):DOUBLE; 

CONST 

NOT.ZERO = IE-30; 

VAR 

A3,C,D,Y_COF^(_COF,LIMIT : DOUBLE; 

BEGIN 

Y_COF := SQR(LONG_AX3S); 
xIcOF := SQR(SHORT_AX3S); 

A := SQR(RHO)*(Y.COF-X_COF); 

B •= 2*X_COF* CENTER. Y*RHO; 

C ; = X cOF*Y_COF-Y_COF*SQR(RHO>X_COF*SQR(CENTER_Y); 

D := (-i+SQRT(SQR(B>4*A*C))/(2*A); 

LIMIT :»ACOS(D); 

IF (LIMIT = PHU) THEN LIMIT := PHI1 + NOT.ZERO; 

PHI.LIMIT := LIMIT; 

END; 


(polynomial extrapolation function ) 

PROCEDURE POLINT(XA,YA : GLN ARRAY; N : INTEGER ;X : DOUBLE ; VAR YJ)Y :DOUBLE); 
VAR 
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NSW : INTEGER; 

W.HP.HO J)IFT,DIFJ)EN : DOUBLE; 
CJ> : GLN ARRAY; 

BEGIN 

NS:=1; 

DEF:«ABS(X-XA[1]); 

FOR I:=l TO N DO 
BEGIN 

DEFT:=ABS(X-XA[I]); 

IF (DIFT < DIF) THEN 
BEGIN 
NS:=I; 

DIF;=DIFT; 

END; 

CfI]:=YA[T]; 

D[I]:=YA[I]; 

END; 

Y:=YA[NS]; 

NS:=NS-1; 

FOR M:=l TO N-l DO 
BEGIN 

FOR I:=l TO N-M DO 
BEGIN 
HO:=XA[I]-X; 

HP:=XA[I+M]-X; 

W:=C[I+1]-D[I]; 

DEN:=HO-HP; 

IF(DEN = 0) THEN 
II := I; 

DEN;=W/DEN; 

D[I]:=HP*DEN; 

C[I]:=HO*DEN; 

END; 

IF (2*NS < N-M)THEN 
DY:=C[NS+1] 

ELSE 

BEGIN 

DY:*>D[NS]; 

NS:=NS-1; 

END; 

Y:=Y+DY; 

END 

END; 


{ integral along phi direction with fixed range in polar coordinate system } 

PROCEDURE TRAPZD_PHI(Y,A,B ; DOUBLE; VAR S : DOUBLE; 

N : INTEGER; VAR GUT : iNTEGER); 

VAR 
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J : INTEGER; 

X,TNM,SUM,DEL : DOUBLE; 


BEGIN 

IF (N = I) THEN 
BEGIN 

S:=0.5*(B-A)*(FUNC(A,Y)+FUNC(B,Y)); 

GUT:=1; 

END 

ELSE 

BEGIN 

TNM:=GLIT; 

DEL:=(B-A)/TNM; 

X:=A+0.5*DEL; 

SUM:=0.0; 

FOR J:=l TO GLIT DO 
BEGIN 

SUM:=SUM+FUNC(X,Y); 

X:-X+DEL; 

END; 

S:=0.5*(S+<B-A)*SUM/rNM); 

GUT:=2+GLIT; 

END 

END; 


( integral along phi direction with fixed range in polar coordinate system } 

FUNCTION QR0MB_PHI(Y,A3 : DOUBLE):DOUBLE; 

LABEL 99; 

CONST 
EPS = 1.0E-6; 

JMAX= 16; 

JMAXP- 17; 

K = 4; 


VAR 

U,GUT_PHI : INTEGER; 

DSS.SS : DOUBLE; 

H,S : ARRAY[1..JMAXP] OF DOUBLE; 
CJ) : GLNARRAY; 

BEGIN 

H[l]:=1.0; 

FOR J;=l TO JMAX DO 
BEGIN 

TRAPZD_PHI(YA3,S[Jy,GUT_PHI); 
IF (J >= K) THEN 
BEGIN 

FOR I := 1 TO K DO 
BEGIN 
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C[I] := H[J-K+I]; 

D[I] := S[J-K+I] 

END; 

POLINT(CJ),K,0.0,SS,DSS); 

IF (ABS(DSS) < (EPS*ABS(SS))) THEN GOTO 99; 
END; 

S[J+1]:=S[J]; 

H[J+1]:=0.25*H[J]; 

END; 

II :»2; 

99: QROMB.PHI := SS; 

END; 


{ integral along range direction on a function whose points are phi_direction 
integrals in polar coordinate system } 

PROCEDURE TRAPZD_RH0(A3 : DOUBLE; VAR S : DOUBLE; 

N : INTEGER; VAR GUT : INTEGER; CIRCLE : INTEGER); 

VAR 

J : INTEGER; 

RHO,TNM,SUM,DELJ»HI2_A,PHI2_B,PHI2_RHO : DOUBLE; 


BEGIN 

IF (N = 1) THEN 
BEGIN 

IF (CIRCLE o 0 ) THEN ( integral on a semi_drcle } 

BEGIN { CIRCLE * 2 for left semi_circle } 

PHI2_A := CIRCLE*PI/2; { CIRCLE = 1 for right semi_cirde ) 
PHI2_B := P HT2 A ( integral limits are PI/2 to PI or 0 to PI/2 } 
END " 

ELSE 

BEGIN 

PHI2_A := Pffl_UMIT(CENTER_Y30NG_AXIS,SH0RT_AXISA); 
PHI2_B := Pffl_LIMIT(CENTER_YJLONG_AXIS,SHORT_AXIS3) 
END; 

R SQRT(SQR(A)+SQR(HEIGHT)); 

R4 := SQR(SQR(R)); 

SIGMA := SIGMA0(ARCTAN(A/HEIGHT)3AR); 
S:=A*QR0MB_PHI(A,PHI1 3HI2_A)*SIGMA/R4; 

R SQRT(SQR(B)+SQR(HEIGHT)); 

R4 := SQR(SQR(R)); 

SIGMA := SIGMA0(ARCTAN(B/HEIGHT)3AR); 

S S+B*QROMB_Pffi(B,PHIU>HI2_B)*SIGMA/R4; 

S:=0.5*(B-A)*S; 

GUT:=1; 

END 

ELSE 

BEGIN 

TNM:=GLIT; 

DEL:*(B-A)/TNM; 
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RHO:= A+0.5 *DEL; 

SUM:=0.0; 

FOR J:=l TO GUT DO 
BEGIN 

IF (CIRCLE o 0 ) THEN 
BEGIN 

Pffl2_RHO :=CIRCLE*PI/2; 

END 

ELSE 

BEGIN ( integral limit on ellipse edge } 

PHI2_RH0 := Pffl_LrMrr(CENTER_Y,LONG_AXIS, SHORT, AXIS.RHO); 
END; 

R := SQRT(SQR(RHO)+SQR(HEIGHT)); { range } 

R4 := SQR(SQR(R)); 

SIGMA := SIGMA 0(ARCTAN(RHO/HEIGHT)J ) AR); { sigmaO | 
SUM:=SUM+RHO*QROMB_PHI(RHO,PHIl J > HI2_RHO)*SIGMA/R4; 
RHO:=RHO+DEL; 

END; 

S:=0.5*(S-HB-A)*SUMTNM); 

GUT:=2*GLIT; 

END 

END; 


PROCEDURE QR0MB_RH0(A3 : DOUBLE; VAR SS : DOUBLE; CIRCLE :INTEGER); 
LABEL 99; 

CONST 
EPS = 1.0E-5; 

JMAX = 16; 

JMAXP = 17; 

K = 4; 

VAR 

U.GLITJRHO : INTEGER; 

DSS : DOUBLE; 

H,S : ARRAY[1.. JMAXP] OF DOUBLE; 

CX > : GLNARRAY; 

BEGIN 

H[l]:=1.0; 

FOR J:=l TO JMAX DO 
BEGIN 

TRAPZD_RHO(A,B,S[Jy,GLIT_RHO, CIRCLE); 

IF (J >= K) THEN 
BEGIN 

FOR I := I TO K DO 
BEGIN 

C[I] := H[J-K+I]; 

DP] :=S[J-K+I] 

END; 
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POLINT(CJ),K,0.0,SS,DSS); 

IF (ABS(DSS) < (EPS+ABS(SS))) THEN GOTO 99; 
END; 

S[J+1]:=S[J]; 

H[J+1]:=0.25*H[J]; 

END; 
n :=3; 

99:; 

END; 


BEGIN 

ASSIGN(FD AT, ’G5_S.DAT’); 

HALFBEAM_INT := 2*BEAMWIDTH_3DB; { integral beam } 

HALFBEAM_EFF := 1.2*BEAMWIDTH_3DB/2; { effective beam : illumination area ) 
S_PAR :* 0.05 *D_R; ( exponential parameter step } 

PAR := 0.8*D_R; ( initial value of exponential para. } 

WHILE (PAR <= (2.5*D_R)) DO 
BEGIN 

APPEND(FDAT); 

WRITELN(FDAT); 

WRITELN(FDATJ > AR/D_R); 

CLOSE(FDAT); 

THETA :* 0; { initial incidence angle } 

WHILE (THETA <51 ) DO 
BEGIN 

II := 0; { II = 0 indicate finishing integral correctly } 

R_THETA := THETA*D_R; 

R0 := HEIGHT/COS(R_THETA); 

Y0 := HEIGHT*TAN(R_THETA); 

ELLIPSE(R_THETA,HALFBEAM_INT,CENTER_YLONG_AXIS,SHORT_AXIS); 
Sl.CIRCLE := 0; 

S2_CIRCLE := 0; 

SINT :* 0; 

CIRCLE := 0; 

PHI1 := 0; 

IF (CENTER_Y < LONG.AXIS) THEN 
BEGIN 

{ integral on left semi_circle } 

CIRCLE := 2; 

PHI I := PI/2; { low integral limit for angle } 

XI :» SHORT_AXIS*SQRT(l-SQR(CENTER_Y/LONG_AXIS)); 

Y1 := LONG_AXIS-CENTER_Y; 

IF (XI > Yl) THEN X_Y :■ XI { up integral limit for range ) 

ELSE X_Y := Yl; 

QROMB_RHO(0.0,X_Y,S2_CIRCLE, CIRCLE); 


( integral on right semi_drcle } 

CIRCLE := 1; 

PHI I := 0; ( low integral limit for angle } 

QROMB_RHO(0.0,X 1 ,S 1_CIRCLE, CIRCLE); 
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CIRCLE := 0; 

RHOl := XI; l low integral limit for range } 

RH02 := CENTER_Y+LONG_AXIS ( up integral limit for range ) 

END 

ELSE 

BEGIN 

RHOl := CENTER_Y-LONG_AXIS; { low integral limit for range 1 
RH02 := CENTER_Y+LONG_AXIS 1 up integral limit for range } 

END; 

IF (RH02 > RHOl) THEN 
{ integral on ellipse } 

QROMB_RHO (RHO 1 ,RH02, SINT, CIRCLE); 

SINT :» SINT+S 1_CIRCLE+S2_CIRCLE; 

fi i .IP SE(R_THETAHALFBEAM_EFF,CENTER_YJ.ONG_AXIS,SHORT_AXIS); 


( calculate illuminated area } 

AILL :» PI *LON G_AXIS *S HORT_AXIS ; 

SINT := 2*SINT*SQR(SQR(R0))/AILL; 

WRrrELN(’PAR,THETA,SIGMAO,SS * ’J»AR:10;5,THETA:10:5,SINT:10:7, 
SINT/SIGMAO(R_THETA r PAR): 10:7); 

IF (II o 0) THEN 

WRTIELN( ’ERROR IN ’ J>AR/D_R:10;5,’ WITH ’.II); 

APPEND(FDAT); 

WRITELN(FDATDB(SINT): 10:6J)B(SINT/SIGMA0<R_THETAJ»AR)); 10:6); 
CLOSE(FDAT); 

FOR JJ := 1 TO 5 DO 
WRITELN; 

THETA := THETA + S.THETA 
END; 

PAR := PAR+S_PAR 
END 
END. 


26 




